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Abstract 



The indirect, charge-carrier mediated coupling between localized magnetic moments is studied for graphene nanoflakes 



of triangular shape and zig-zag edge. The characteristic feature of such nanoflakes is the presence of a shell of zero- 
energy states in the electronic spectrum. The tight-binding Hamiltonian supplemented with a Hubbard term is used 

t^- for electronic structure calculations. The indirect RKKY (Ruderman-Kittel-Kasuya-Yosida) coupling energy is de- 

rived from the total electronic energy of the system in a non-perturbative way. The attention is focused on the on-site 
and plaquette impurities situated along the edge. The charge doping is also taken into account. It is found that the 
zero-energy states may give rise to a coupling mechanism which is describable by a first-order perturbation calculus 
and can yield robust indirect coupling of both ferro- and antiferromagnetic character, which dominates ovet the usual 

i-G RKKY mechanism. The numerical results obtained emphasize the importance of the Hubbard term and the effect of 

r/} the charge doping on the coupling. 
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1. Introduction 

Graphene, a first two-dimensional material 01121' opens the route to new physics and applications. Among many 
others, prospects for graphene-based spintronics [3| encourage extensive studies of magnetic properties of graphene 
and its derivatives @). As a consequence, graphene nanostructures focus considerable efforts of researchers |5|. 
Graphene nanoflakes behave in a quantum-dot manner, exhibiting discretization of energetic spectrum [6|, what mod- 



pvi ifies their behaviour with respect to those expected for infinite graphene monolayers and allows for emergence of new 

phenomena. In particular, graphene nanoflakes of triangular shape attracted particular attention [7-231, mainly due to 

t— I their peculiar electronic spectrum, containing a shell of degenerate zero-energy states leading to nontrivial magnetic 

properties for zigzag-edged structures. 

■^J" One of the issues related to magnetism in graphene systems is the problem of an indirect coupling between 

localized magnetic moments in graphene, mediated by the charge carriers [24-47|. This kind of coupling, known 
as Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction l48W50ll . presents an unique behaviour in graphene, due to 
its linear dispersion relation, unlike this found in metallic two-dimensional systems [51]. Moreover, the oscillatory 
properties of RKKY coupling between on-site impurities in graphene are ruled by the principle relating the sign of 
coupling integral to the sublattice in which both magnetic impurities are placed. On the other hand, in geometrically 
confined graphene systems, the indirect coupling differs from that in infinite systems B31l52l454l . The already studied 
systems include nanoflakes with zig-zag and armchair edges and equal number of carbon atoms in two sublattices 
ll52l as well as graphene nanoribbons of both kinds of edges ll24l l45l 1531 . Let us mention that in such circumstances, 
some additional mechanisms of indirect exchange, dissimilar to the typical RKKY mechanism, may become operative 
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M =3 hexagonal rings 

Figure 1 : Schematic view of a triangular zigzag-edged nanoflake with M = 3 hexagonal rings at each side; various impurity pairs considered in the 
calculations are marked; empty (filled) circles denote atoms belonging to the majority (minority) sublattice. 

Il52"ll53l . This fact, together with the unique electronic properties of zigzag triangular graphene nanoflakes, makes 
them particularly interesting systems for further studies of indirect coupling. 

Moreover, it can be mentioned that the magnetically doped nanosized objects provide the opportunity of control 
over the magnetic properties by means of doping with a small number of the charge carriers . This involves both diluted 
magnetic semiconductor quantum dots 11551 1561 and graphene-based nanostructures ll52l and provides an additional 
reason for studies of such nanoscopic systems. 

The aim of the present paper is to discuss the form of an indirect coupling between the magnetic impurities in 
graphene nanoflakes of triangular shape and zigzag edge. The calculations are based on the tight-binding approxima- 
tion with Hubbard term for the electronic structure. The charge doping of the nanostructures is taken into account. 

2. Theory 

The subject of the study are the triangular graphene nanoflakes posessing zig-zag edge. The size of such a flake 
can be described using the number M of hexagonal rings forming each side of the structure (see Fig. fTh. Such GNF 
consists of N = M 2 +4M + 1 carbon atoms [8|, and the two sublattices contain unequal number of atoms. Therefore, 
one sublattice (empty circles) is the majority sublattice, while the other sublattice (filled circles) is the minority one. 
The difference on the number of atoms between them amounts to M - 1 . 

In the present study, the electronic structure of the GNFs is described using the following Hamiltonian for p z 
electrons of carbon atoms [52]: 



('J>.0~ ' 



The summation over (i, j) denotes the nearest-neighbour pairs of sites, while cr =1, f is the electron spin. The total 
Hamiltonian involves a tight-binding nearest-neighbours hopping term with hopping integral t (approximately equal to 
2.8 eV 11571 ). In addition, the coulombic interactions between electrons are described by means of an on-site Hubbard 
term, to which a mean field approximation (MFA) is applied in a self-consistent way J4] [25] [52] [58] [59] . Let us 
mention that, according to [59], the MFA results for total energy compare successfully to exact studies of the Hubbard 
Hamiltonian for graphene nanostructures provided that approximately U/t < 2, what justifies its application here. 

It is worth mentioning that the value of U used in such calculations does not represent an actual on-site Coulombic 
energy. Instead, it is an effective parameter influenced by the long-range nature of Coulombic interactions and its value 
for small graphene nanosy stems should be close to U/t ^ 1 [60]. 

It follows from the studies of the electronic structure (see for example JSjITOl or lfTTl[T3l ) that a zigzag nanoflake 
with M hexagons at the edge possesses a shell of M - 1 (degenerate) zero-energy states in its spectrum. The most 
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Figure 2: (a) Total ground-state spin of the nanoflakes with M = 2, 3,4 as a function of charge-carrier doping; (b) total spin associated with the 
lattice sites of the nanoflake with M = 3 for charge neutrality in the ground-state, for U/t — * and U/t = 1; (c) Absolute value of the spin density 
associated with various lattice sites of the M = 3 nanoflake as a function of the normalized Hubbard on-site exchange energy U/t. 
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typical situation is that of a charge neutral, undoped nanoflake, with Aq — (when the number of the electrons in the 
flake is equal to the number of carbon atoms, what corresponds to the state of half-filling of the energy levels). In such 
a situation the zero-energy states lie at the Fermi level. Let us mention that in presence of the coulombic interaction 
between the electrons, it is energetically favourable to fill these states according to the Hund rule, so that at half-filling 
all the electrons occupying the zero-energy states have parallel spins. Therefore, the state with Aq - leads to a 
maximum spin of the nanoflake, equal to (M— l)/2 Bl [TTll6Tll62l . Doping with either electrons (i.e. Aq > 0) or holes 
{Aq < 0) leads to decrease of the net magnetic moment, which vanishes for completely empty or completely filled 
shell of zero-energy states. Once the shell is either empty or completely filled, further doping may cause the magnetic 
moment to increase again. The dependence of the total spin of the nanoflake on the charge doping resulting from our 
Hamiltonian is illustrated for three nanoflakes, M = 2, 3,4, in Fig.pja). Let us note that the doped triangular flakes 
were subject to extensive studies in ATI [TBI l63l . 

The existence of the zero-energy states and the form of their wavefunctions can be expected to have a crucial 
impact on the indirect coupling between the magnetic moments in the nanoflakes. In particular, the nonzero spin 
distribution arising due to the presence of zero-energy states influences the mentioned properties. 

Let us illustrate briefly the spin distribution associated with the zero-energy states of zigzag triangular nanoflakes 
(for an example of the individual states in a nanoflake with M — 4 see Fig. 5 in [ 10]). In Fig.|2tb), the total spin density 
(which corresponds to the difference between the charge density of spin up and spin down electrons) for the nanoflake 
with M — 3 is shown for Aq — and in the limit of U/t — > 0. The radii of the circles correspond proportionally 
to the spin magnitude at the given lattice site while different colours depict the opposite signs of spin density. It is 
visible that the nonzero spin density is present only at the majority sublattice sites. Moreover, the wavefunctions of 
the zero-energy states are strongly localized at the zigzag edges of the nanoflake. The situation is changed when finite 
U/t is introduced. For U/t - 1 [Fig[2jb)], the spin density takes the opposite signs at both sublattices, however, it is 
still dominant at the minority sublattice sites. The spin orientations at both sublattices are opposite. The dependence 
of the spin density at four representative lattice sites at the zigzag edge is plotted in Fig. |2|c) for a wide range of the 
Hubbard exchange energies U/t. It is visible that the presence of coulombic interactions strongly promotes the spin 
polarization of both sublattices (of course conserving the total nanoflake spin). The increase rate of the polarization 
tends to drop for large U/t, reaching the value of a few, when the spin associated with each site reaches noticeably 
high values. Let us mention that a study of the magnetization of the finite graphene nanoflake as a function of the U/t 
has been performed in |59| and a similar increase has been found within mean-field results. Comparison with exact 
diagonalization-based calculations shows that the rise of magnetization with Hubbard exchange energy is rather more 
linear-like. However, the results for U/t < 2 obtained within mean-field approach are quite accurate (this range located 
the system below a critical value of U/t for transition to Mott-Hubbard insulator). Let us explain also that there is no 
spontaneous spin polarization of the nanoflake strictly for U/t = 0, unlike it happens in the limit of U/t — > 0, when 
the vanishing coulombic interactions still break the symmetry [Fig. |2jb)]. However, when the symmetry is broken 
for example by the presence of a pair of magnetic impurities, the results for weak U/t are convergent to the system 
behaviour for U/t = 0. Therefore, in further presentation of the indirect coupling calculations distinguishing between 
the cases of U/t = and U/t — > is not critical. 

Then let us consider a nanoflake with a pair of localized magnetic impurity moments. In presence of these mo- 
ments, the total Hamiltonian 'H of the system takes the form of 'Ho + < H mp . Here, < Hi mp is the Anderson-Kondo term 
for the interaction between the magnetic impurity spins and spins of the charge carriers, with exchange energy equal 
to /. For on-site impurities located at the sites a, b, the term is of the form < Hi, np - iS\ {n a ^ - n a \) + ^S z h {rif,^ - ra^j,), 
while for plaquette impurities, positioned in the centers of the hexagons, it yields < Hi mp = lS z a iZieNN i( n i,t ~ n Ui)\ + 
5^i YiieNN i( n i,1 ~ n u)] ( wnere the summation is over all the sites belonging to the hexagon housing the impurity). 

The total electronic energy E, ot is then obtained for a nanoflake with N+Aq electrons, separately for ferromagnetic 
(F) and antiferromagnetic (AF) orientation of the magnetic impurity spins. This allows for determination of an indirect 
exchange coupling RKKY-like integral from 2S 2 J RKKY = Ej£ - Ef ot E4ll52l . Therefore, the interaction between 
the impurity spins is projected onto an effective Hamiltonian of the form H a t, = -J RKKY S z a S z b . The Ising form 
of Anderson-Kondo term leads to Ising-like interaction between impurity spins, however, for isotropic Heisenberg 
exchange the Hamiltonian H a i, would be of the Heisenberg form with coupling energy J RKKY unchanged. Therefore 
the choice of the Ising form is only for the sake of facilitating the numerical calculations. Let us note that, due to 
preserved electron-hole symmetry of the total Hamiltonian, all the results are expected to be independent on the sign 
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Figure 3: Indirect coupling energy values as a function of charge-carrier doping for two selected nanoflakes, M = 2, 4, in absence (U/t = 0) and 
in presence (U/t = 1) of the Hubbard term, (a) 2A 2nd neighbours pair; (b) 2B 2nd neighbours pair; (c) 3rd neighbours pair; (d) 1st neighbours 
plaquette pair. For explanation see Fig.fT] 

of Aq, i.e. identical for electron- and hole-doping. We also must emphasize that the present method yields the total 

energy of indirect coupling, which may involve mechanisms which differ from the usual RKKY interaction mechanism 

operative via electronic states occupied by pairs of opposite spin electrons. However, let us call the resulting coupling 
jRkky for b rev i t y 

3. Numerical results 



The numerical results, presented in the subsequent figures were obtained by means of the total energy calculation, 
on the basis of exact diagonalization of the Hamiltonian (fTl in the single particle approximation. For numerical 
calculations, the LAPACK package was used [64]. Let us note that these results are non-perturbational with respect to 
the J parameter, describing the interaction between impurity spins and spins of the charge carriers, like in [ 24 . 22 [52) . 
This allows to capture the effects beyond the second-order perturbational calculus, the presence of which can be 
expected in such molecular- like systems; on the other hand, only the numerical results can be obtained. For the 
further studies, the value of J/t = 0.1 was accepted (unless the /-dependence of coupling energy is studied). 

The triangular GNFs considered in the present study were characterized by the values of M ranging between 
M — 2 and 7. Several representative placements of the magnetic impurities were taken into consideration, strongly 
focusing on the position close to the zig-zag edge. For on-site impurities, two inequivalent second-neighbour pairs, 
2A and 2B (see Fig. fTb were discussed (i.e. impurities in minority or in majority sublattice) in as well as the third- 
neighbour impurities placed in different sublattices. In addition, plaquette impurities (placed in the centres of the 
hexagons) were also considered in nearest-neighbour position. For studies of the distance dependence of the indirect 
coupling, also the other impurity positions were taken into account [as illustrated in Fig. 171 a)]. 

At the beginning, the importance of the Hubbard term for the indirect exchange integrals is verified. Fig.[3]presents 
a comparison of the values of coupling energies obtained for U/t = and for U/t = 1, the latter being the value 
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Figure 4: Indirect coupling energy values as a function of the localized spin-electron spin exchange energy /, for charge neutrality Aq = and for 
two values of doping \Aq\ = 1,2, for a nanoflake with M = 3. (a) 2A 2nd neighbours pair; (b) 2B 2nd neighbours pair; (c) 3rd neighbours pair; (d) 
1st neighbours plaquette pair. For explanation see Fig.fT] 



accepted for the Hubbard term for further calculations. The dependence of the indirect coupling on the charge doping 
Aq is demonstrated for two flakes characterized by M — 2 and M - 4. The results are plotted for 2A and 2B second- 
neighbour on-site impurities, for third-neighbour on-site impurities and for nearest-neighbour plaquette impurities, 
all located at the edge of the nanoflake. In general, it can be observed that for vanishing coulombic interactions, the 
indirect coupling is very weak for the cases of 2nd neighbours in minority sublattice or for 3rd neighbours, close to the 
charge neutrality. It remains negligible as long as the zero-energy states shell is partially filled (i.e., for \Aq\ < M - 1). 
On the contrary, for plaquette impurities and for 2nd neighbours in majority sublattice, the interaction has considerable 
magnitude in the vicinity of the charge neutrality point, for \Aq\ < M - 1 . Including the Hubbard term with U/t - 1 
crucially influences the behaviour of an indirect coupling for 2A and 3rd neighbours position. Namely, the interaction 
energy is increased (toward ferromagnetic values for 2nd neighbours and antiferromagnetic values for 3rd neighbours). 
This effect remains robust against charge doping until \Aq\ < M - 1. For 2B and plaquette impurities, the effect is 
similar in its direction, however, much less pronounced. In all cases, the peak coupling energies correspond to Aq = 0. 
Let us observe another peaks for indirect coupling present at \Aq\ = 4 for M - 2, corresponding to another maximum 
nanoflake spin configuration [see Fig.[2|a)]. 

In order to explain this behaviour, let us focus on the dependence of the indirect coupling energies on the exchange 
parameter J, describing the interaction between localized spin and electron spin. Such results are illustrated in Fig. HI 
for M - 3 and four impurity locations at the nanoflake edge, in doubly logarithmic scale. Our calculations rely on 
non-perturbative procedure of determining the total energy, thus are able to identify the contributions to an indirect 
coupling described by various orders of perturbational calculus. For U/t = and for 2A and 3rd neighbours, the 
coupling varies quadratically with J and is weak in magnitude. This suggests the origin from the second-order 
perturbational process, which is the typical origin of indirect RKKY coupling. In the mentioned cases, the interaction 
is also weakly dependent on the charge doping. If U/t - 1 is taken into consideration, the situation is changed. Then, 
the indirect coupling becomes a linear function of |/| and the magnitude is greatly pronounced, if \Aq\ < 2, while 
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Figure 5: Indirect coupling energy values as a function of charge-carrier doping for increasing size of the nanoflake M = 2, 3, 5, 7, in presence of 
the Hubbard term (U/t = 1). (a) 2A 2nd neighbours pair; (b) 2B 2nd neighbours pair; (c) 3rd neighbours pair; (d) 1st neighbours plaquette pair. 
For explanation see Fig.fT] 



its behaviour resembles the case of U/t - for \Aq\ - 2. Such a behaviour indicates, that a leading contribution 
originates form the first-order perturbational process, involving the zero-energy states, the shell of which is partially 
filled for \Aq\ < 2. On the other hand, for 2B neighbours and plaquette impurities, both for U/t = and 1 the linear 
dependence of indirect coupling on |/| is observed for |A<gr| < 2, and again, for empty shell of the zero-energy states the 
usual quadratic behaviour is restored. As it can be deduced form the Fig. |2jb) and (c), nonzero spin density coming 
from a half-filled shell of the zero-energy states is present only at the majority sublattice for U/t — > 0, and switching 
on finite U/t causes the redistribution of the spin density to both sublattices. Therefore, for U/t — > 0, the first-order 
perturbational contribution is possible only for the impurities in the majority sublattice (or for plaquette impurities, 
each of which interacts with equal number of electrons residing at both sublattices). Increase of U/t enables the 
same mechanism for impurity pairs in minority sublattice (e.g. 2A) or for impurities in both sublattices (e.g. 3rd 
neighbours). 

Let us consider the first-order perturbational mechanism quantitatively |52|. The first-order correction to energy 
of all the states due to presence of impurity spin S a in site a is AE = -\J\S z a s z , and is proportional to the total spin in 
this site s z a = (n a j - n a ^) /2. Therefore, if the spin-degenerate states are filled with pairs of opposite-spin electrons 
and no net spin polarization is present, the first order correction vanishes. Let us assume that two impurity spins 
are present at the sites a and b (assuming \s z a \ > \s\\). If the impurity spins are polarized ferromagnetically, the total 
energy change is AE^ b - -\J\S s z a - t]\J\S s z b , where rj - s z a s z b / \s z a s z b \ is the sign of s z a s z b . The orientation of the spins s z a 
and s z b depends on the sign of J such that it maximizes the absolute value of the correction (which itself is negative), 
of course conserving the relative orientation of s z a and s z b . On the other hand, for antiferromagnetic orientation of 
impurity spins, it yields AEf ib = -\J\S s z a + rj\J\S s z h . The indirect coupling energy due to this mechanism is then 
2S 2 J RKKY - AE^ b - AE^ b = 2rjS\J\s z b . Let us note that owing to linear dependence on \J\, it dominates over the usual 
J 2 -proportional mechanism. 
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The idea of the first-order perturbational mechanism can be expressed briefly by the statement that the energetically 
favourable mutual orientation of impurity spins follows the mutual orientation of the electron spin at sites a and b. 
Therefore, it is ferromagnetic for impurities in the same sublattice and antiferromagnetic in different sublattices [see 
Fig. gb)]. 

In Fig. Hit is visible, that for strong exchange coupling between impurity spins and electron spins J/t, taking the 
value of a few, the coupling energy tends to saturate, and then even decrease; also the change of sign is possible in 
some cases, as for 2B 2nd neighbours and plaquette impurities, in presence of Hubbard term and for empty zero- 
energy shell. 

The evolution of an indirect coupling with the size of triangular nanoflake can be analysed in Fig. B] where the 
dependence of coupling energy on charge doping is presented for M - 2, ... ,7. The results take into account the 
value of U/t — 1. In all the presented cases, it is visible that the coupling depends in principle rather weakly on 
charge doping as long as \Aq\ < M - 1, i.e. as long as the shell of zero-energy states is neither completely empty 
nor completely filled. The coupling energy in general gradually decreases with doping within mentioned range. The 
situation is slightly different only for M - 7, where for 2A and 3rd neighbours, the coupling diminishes at \Aq\ = 3. 
It can be concluded that a quite strong indirect coupling, robust against charge doping, can be achieved at the zigzag 
edge of the triangular nanoflake. 

The effect of the Hubbard exchange energy parameter U on the indirect coupling energy can be followed in 
the Fig. [6| where the cases of two 2nd neighbour impurities as well as 3rd neighbours and plaquette impurities are 
studied. The state of the charge neutrality and moderate dopings (\Aq\ = 1,2) are considered. For 2nd neighbours 
2A [Fig. |6ja)] and for 3rd neighbours [Fig.lfjc)], the coupling energy increases noticeably with the increase of U, 
strengthening the interaction over orders of magnitude. This is due to the mentioned fact, that for U/t - 0, there is 
no first-order contribution to indirect coupling for these impurity pairs (as the spin density vanishes at the minority 
sublattice under such circumstances). As Hubbard exchange is switched on, the opposite-sign spin density develops 
at both sublattices [see Fig.|2jb) and (c)] and first-order contribution (of ferro- or antiferromagnetic kind) appears and 
increases in strength. When a zero-energy states shell is empty (\Aq\ - 2), the changes in indirect coupling energy are 
much less pronounced. The situation is somehow different for 2B 2nd neighbours and plaquette impurities [Fig.|6tb) 
and (d), respectively], where first-order contribution to indirect coupling exists already at U/t - 0, so that the increase 
in Hubbard exchange energy only moderately influences the interaction between impurities for \Aq\ < 2. For empty 
zero-energy states shell, the indirect coupling tends to change sign when U/t increases, which can be attributed to 
charge density redistribution. This happens also under the influence of increasing J/t [see Fig.[4|b)]. Note that the 
coupling magnitude in these cases is greatly reduced in comparison with 2 A 2nd neighbours or 3rd neighbours. 

It should be emphasized that for plaquette impurities in undoped graphene monolayer, their indirect interactions 
are expected to be antiferromagnetic ll24l . However, at the edge of the studied nanoflakes, the coupling has ferro- 
magnetic sign for half-filling conditions and for low charge dopings, as shown in Fig.[6|d). This is mainly due to the 
influence of the first-order coupling mechanism. The spin densities originating from zero-energy states are opposite 
in sign, but of unequal absolute values for both sublattices [see Fig. Efb)], so that a nonzero net spin polarization is 
present at a given hexagon, being especially pronounced for locations close to the edge. Moreover, the orientation of 
net spin for each hexagon is the same. Therefore, a ferromagnetic coupling is established between plaquette impurities 
under such conditions. This constitutes an important finite-size effect in the indirect coupling in the studied structures. 

Finally, let us study the distance dependence of an indirect coupling along the zigzag edge of an undoped 
nanoflake. First we focus on the on-site magnetic impurities and for the purpose of this analysis, one of the spins 
is located in majority sublattice site closest to the edge, while the other takes various positions along the zigzag edge, 
starting from nearest neighbours [as illustrated in Fig. [7Ta)]. The calculations are performed for a nanoflake with 
M -1 . The results are depicted as a function of distance in Fig. |7|b) for three values of U/t; 0, 0. 1 and 1 . It is visible 
that without coulombic interactions included, the coupling tends to vanish considerably fast (it can be verified that the 
decay is approximately exponential in character or even faster). Its sign oscillates according to the rule that coupling 
is ferromagnetic for on-site impurities in the same sublattice while it is antiferromagnetic for different sublattices for 
half-filling conditions. Inclusion of the Hubbard term significantly changes the distance behaviour of an indirect cou- 
pling. For U/t = 0.1, the coupling magnitudes are enhanced and the distance dependence of coupling becomes rather 
weak after some initial fall. For U/t - 1, the coupling magnitude becomes almost distance-independent (regarded 
separately for the same and different sublattices) and the sign oscillates between the two cases. This reflects the 
dominance of the discussed coupling mechanism based on zero-energy states. For significant U/t value, the coupling 
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Figure 6: Indirect coupling energy values as a function of the Hubbard on-site exchange energy U/t, for charge neutrality Aq = and for two 
values of doping \Aq\ = 1,2, for a nanoflake with M = 3. (a) 2A 2nd neighbours pair; (b) 2B 2nd neighbours pair; (c) 3rd neighbours pair; (d) 1st 
neighbours plaquette pair. For explanation see Fig.fT] 

energies are weakly dependent on the distance owing to almost uniform distribution of spin density values along the 
edge for both sublattices [see Fig.[2tb)]. 

Secondly, we study the distance dependence for the case of plaquette impurities located along the edge [see 
Fig.[7|a)]. The results are presented in Fig.[7|c) in logarithmic scale. The coupling is always ferromagnetic along 
the edge and in absence of Hubbard term it decays exponentially with the distance for small distances, while the 
decay becomes even faster for larger separation of impurities. Switching on the coulombic interactions enhances 
the coupling energy and flattens the distance dependence. Like for on-site impurities, for U/t - 1, the coupling 
loses its distance dependence along the zigzag edge. Once more, let us mention that in this case the coupling sign is 
ferromagnetic, while the opposite sign is predicted for infinite graphene monolayer. 

Let us note that an influence on the coulombic interactions taken into account by means of a Hubbard term has been 
studied for monolayer graphene and zigzag nanoribbons by Black-Schaffer [25]. In particular, a distance-independent 
coupling at zigzag edge of a narrow zigazg nanoribbon (which indicated spontaneous magnetic polarization at the 
edge) has been found in accordance with the present findings. Moreover, also the mechanism of a first-order con- 
tribution to the coupling, proportional to the electron spin at one of the impurity sites has been identified in such 
case. 



4. Concluding remarks 

In the paper, the indirect coupling between localized on-site and plaquette magnetic moments on the triangular 
graphene nanoflakes with zig-zag edge was studied. Particular attention was paid to the role of a shell of zero-energy 
states present in these nanostructures. In that context, the influence of charge coupling has been studied and the 
importance of Hubbard on-site exchange term was emphasized. In the presence of edge-localized zero-energy states, 
the placement of the impurities at the nanoflake edge yields particularly interesting behaviour of an indirect coupling. 
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(c) Indirect coupling energy values as a function of the distance between the plaquette impurities for a nanoflake with M = 1 for various strengths 
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Let us note that the edge-localized states for example in pits in graphene layer have been observed by means of 
STM/STS methods already in [65, 66 1. Due to the edge-localized form of the wavefunction of the zero-energy states 
and nonzero spin density at the ground state for charge neutrality, a first-order perturbational mechanism becomes 
operative for many impurity positions along the edge. This kind of an indirect coupling, similar to a double exchange, 
dominates strongly over the usual RKKY second-order perturbative mechanism in numerous cases. Depending on 
the on-site impurity location, the enhanced indirect coupling can be either ferromagnetic (for impurities in the same 
sublattice) or antiferromagnetic (impurities in various sublattices). If the exchange energy between localized on-site 
spins and electron spins is not too large, we find no deviation from this sign rule related to RKKY coupling in graphene 
for on-site impurities. The situation is different for plaquette impurities. While in undoped graphene monolayer the 
RKKY coupling between such impurities is antiferromagnetic [241 . our results show that it becomes ferromagnetic 
for impurities located at the edge of a triangular nanoflake. This sign change is attributed mainly to the first-order 
coupling mechanism. 

The partially filled shell of zero-energy states guarantees the coupling which is quite robust against charge doping 
in larger nanoflakes; on the other hand, in the smallest structures it can be switched form large value to a small one by 
doping with a single charge carrier only. 

The distance dependence of the coupling along the edge is also heavily influenced by the presence of the coulom- 
bic on-site interactions. Namely, the fast decay of the coupling for U/t - is replaced with almost flat distance 
dependence (distinguishing only between the sublattices) for U/t = 1 . 

The obtained results for triangular nanoflakes demonstrate that the indirect coupling between impurity spins in 
such system shows a nontrivial dependence on the charge doping and impurity position. This properties originate 
mainly from a subtle interplay between first- and second-order perturbative contributions. Such a rich behaviour may 
motivate further studies of routes to control the coupling. Moreover, it might also stimulate further works on spin 
readout in such molecular magnet-like systems (see e.g. Il67l ). 

5. Acknowledgments 

This work has been supported by Polish Ministry of Science and Higher Education on a special purpose grant 
to fund the research and development activities and tasks associated with them, serving the development of young 
scientists and doctoral students. 

The computational support on HUGO cluster at Department of Theoretical Physics and Astrophysics, P. J. Safarik 
University in Kosice is gratefully acknowledged. 

The author is deeply grateful to T. Balcerzak for critical reading of the manuscript. 

References 
References 

[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva. A. A. Firsov, Science 306 (2004) 666. 

[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov, Nature 438 (2005) 
197. 

[3] C. Jozsa, B. J. van Wees, Graphene spintronics, in: E. Y. Tsymbal, I. Zutic (Eds.), Handbook of Spin Transport and Magnetism, Chapman 
and Hall/CRC, 201 1, p. 579. 

[4] O. V. Yazyev, Reports on Progress in Physics 73 (2010) 056501. 

[5] I. Snook, A. Barnard, in: S. Mikhailov (Ed.), Physics and Applications of Graphene - Theory, InTech, 201 1, Ch. 13, p. 277. 

[6] M. Wiessner, N. S. R. Lastra, J. Ziroff. F. Forster, P. Puschnig, L. Dossel, K. Mllen, A. Scholl, F. Reinert, New Journal of Physics 14 (2012) 
113008. 

[7] M. Ezawa, AIP Conference Proceedings 1399 (2011) 775. 

[8] M. Ezawa, Journal of Nanoscience and Nanotechnology 12 (2012) 386. 

[9] M. Ezawa, Physical Review B 76 (2007) 245415. 
[10] M. Ezawa, in: S. Mikhailov (Ed.), Physics and Applications of Graphene - Theory, InTech, 201 1, Ch. 12, p. 243. 
[11] P. Potasz, A. D. Giiclii, A. Wojs, P. Hawrylak, Physical Review B 85 (2012) 07543 1 . 
[12] O. Voznyy, A. D. Giiclii, P. Potasz, P. Hawrylak, Physical Review B 83 (201 1) 165417. 
[13] P. Potasz. A. D. Giiclii, P. Hawrylak, Physical Review B 81 (2010) 033403. 
[14] A. Zhou, W. Sheng, Journal of Applied Physics 112 (2012) 094313. 
[15] W. Sheng, M. Korkusinski, A. D. Giiclii, M. Zielinski, P. Potasz, E. Kadantsev, O. Voznyy, P. Hawrylak, Frontiers of Physics 7 (2012) 328. 

11 



NOTICE: this is the authors version of a work that was accepted for publication in Physica E: Low-dimensional 
Systems and Nanostructures. Changes resulting from the publishing process, such as peer review, editing, 
corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. 
Changes may have been made to this work since it was submitted for publication. A definitive version was 
subsequently published in Physica E: Low-dimensional Systems and Nanostructures, [VOL 52, (August 2013)] 
DOI:10.1016/j.physe.2013.03.017 

[16] X. Li, J. Zhou, Q. Wang, X. Chen, Y. Kawazoe, P. Jena, New Journal of Physics 14 (2012) 033043. 

[17] H. P. Heiskanen, M. Manninen, J. Akola, New Journal of Physics 10 (2008) 103015. 

[18] J. Fernandez-Rossier, J. J. Palacios, Physical Review Letters 99 (2007) 177204. 

[19] M. Zarenia, A. Chaves, G. A. Farias, F. M. Peeters, Physical Review B 84 (201 1) 245403. 

[20] S. Zhang, J. Zhou, X. Li, Q. Wang, Journal of Nanoparticle Research 14 (2012) 1. 

[21] M. R. Philpott, F. Cimpoesu, Y. Kawazoe, Chemical Physics 354 (2008) 1. 

[22] B. Uchoa, T. G. Rappoport, A. H. Castro Neto, Physical Review Letters 106 (2011) 016801. 

[23] A. M. Silva, M. S. Pires, V. N. Freire, E. L. Albuquerque, D. L. Azevedo, E. W. S. Caetano, The Journal of Physical Chemistry C 1 14 (2010) 

17472. 

[24] A. M. Black-Schaffer, Physical Review B 81 (2010) 205416. 

[25] A. M. Black-Schaffer, Physical Review B 82 (2010) 073409. 

[26] V. K. Dugaev, V. I. Litvinov, J. Barnas, Physical Review B 74 (2006) 224438. 

[27] L. Jiang, X. L, W. Gao, G. Yu, Z. Liu, Y Zheng, Journal of Physics: Condensed Matter 24 (2012) 206003. 

[28] F. Peng, W. Hongbin, Physica B: Condensed Matter 407 (2012) 34343436. 

[29] S. R. Power, F. S. M. Guimaraes, A. T. Costa, R. B. Muniz, M. S. Ferreira, Physical Review B 85 (2012) 19541 1 . 

[30] M. Sherafati, S. Satpathy, Physical Review B 83 (2011) 165425. 

[31] S. R. Power, M. S. Ferreira, Physical Review B 83 (2011) 155432. 

[32] H. Lee, J. Kim, E. R. Mucciolo, G. Bouzerar, S. Kettemann, Physical Review B 85 (2012) 075420. 

[33] H. Lee, E. R. Mucciolo, G. Bouzerar, S. Kettemann, International Journal of Modern Physics Conference Series 1 1 (2012) 177. 

[34] S. R. Power, P. D. Gorman, J. M. Duffy, M. S. Ferreira, Physical Review B 86 (2012) 195423. 

[35] S. R. Power, M. S. Ferreira, Crystals 3 (2013) 49. 

[36] M. Sherafati, S. Satpathy, Physical Review B 84 (201 1) 125416. 

[37] E. Kogan, Physical Review B 84 (201 1) 1 15 1 19. 

[38] E. Kogan, Graphene 2 (2013) 8. 

[39] L. Brey, H. A. Fertig, S. Das Sarma, Physical Review Letters 99 (2007) 1 16802. 

[40] M. Killi, D. Heidarian, A. Paramekanti, New Journal of Physics 13 (2011) 053043. 

[41] S. Saremi, Physical Review B 76 (2007) 184430. 

[42] J. E. Bunder, H.-H. Lin, Physical Review B 80 (2009) 153414. 

[43] V. V. Cheianov, The European Physical Journal - Special Topics 148 (2007) 55. 

[44] B. Wunsch, T Stauber, F. Sols, F. Guinea, New Journal of Physics 8 (2006) 318. 

[45] J. Klinovaja, D. Loss, Physical Review B 87 (2013) 045422. 

[46] F. Parhizgar, R. Asgari, S. H. Abedinpour, M. Zareyan, Physical Rreview B 87 (2013) 125402. 

[47] O. Roslyak, G. Gumbs, D. Huang, Journal of Applied Physics 1 13 (2013) 123702. 

[48] M. A. Ruderman, C. Kittel, Physical Review 96 (1954) 99. 

[49] T. Kasuya, Progress of Theoretical Physics 16 (1956) 45. 

[50] K. Yosida, Physical Review 106 (1957) 893. 

[51] M. T Beal-Monod, Physical Review B 36 (1987) 8835. 

[52] K. Szatowski, Physical Review B 84 (2011) 205409. 

[53] K. Szalowski, Journal of Physics: Condensed Matter 25 (2013) 166001. 

[54] J. E. Bunder, J. M. Hill, The Journal of Chemical Physics 136 (2012) 154504. 

[55] J. Fernandez-Rossier, L. Brey, Physical Review Letters 93 (2004) 117201. 

[56] F. Qu, P. Hawrylak, Physical Review Letters 96 (2006) 157201. 

[57] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim, Reviews of Modern Physics 81 (2009) 109. 

[58] M. Fujita, K. Wakabayashi, K. Nakada, K. Kusakabe, Journal of the Physical Society of Japan 65 (1996) 1920. 

[59] H. Feldner, Z. Y. Meng, A. Honecker, D. Cabra, S. Wessel, F. F. Assaad, Physical Review B 81 (2010) 1 15416. 

[60] M. Schiiler, M. Rosner, T O. Wehling, A. I. Lichtenstein, M. I. Katsnelson, |arXiv:1 302 . 1437 

[61] W. L. Wang, S. Meng, E. Kaxiras, Nano Letters 8 (2008) 241. 

[62] M. Ezawa, Physica E: Low-dimensional Systems and Nanostructures 40 (2008) 1421. 

[63] A. D. Giiclu, P. Potasz, O. Voznyy, M. Korkusinski, P. Hawrylak, Physical Review Letters 103 (2009) 246805. 

[64] E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, 

D. Sorensen, LAPACK Users' Guide, 3rd Edition, Society for Industrial and Applied Mathematics, Philadelphia, PA, 1999. 

[65] Z. Klusek, Vacuum 63 (2001) 139. 

[66] Z. Klusek, W. Kozlowski, Z. Waqar, S. Datta, J. Burnell-Gray, I. Makarenko, N. Gall, E. Rutkov, A. Tontegode, A. Titkov, Applied Surface 

Science 252 (2005) 1221. 

[67] A. Candini, S. Klyatskaya, M. Ruben, W. Wernsdorfer, M. Affronte, Nano Letters 1 1 (201 1) 2634. 



12 



